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Abstract: Calculations for fission and cluster decay of Sr are presented for this nucleus 
to be in its ground-state or formed as an excited compound system in heavy-ion reactions. 
The predicted mass distribution, for the dynamical collective mass transfer process assumed 
for fission of 76 SV, is clearly asymmetric, favouring a-nuclei. Cluster decay is studied within 
a preformed cluster model, both for ground-state to ground-state decays and from excited 
compound system to the ground-state(s) or excited states(s) of the fragments. 

JJSV is a superdeformed nucleus with an estimated quadrupole deformation 02 = 0.44 (see 
Fig. 1 in [1]). From the point of view of known spherical shell closures at Z=N=40, such a 
large ground-state deformation for 76 Sr means the natural breaking of these spherical shells 
and hence nuclear instability against both fission and exotic cluster decay processes. However, 
we shall see in the following that, like the other superdeformed nuclei in this mass region [1], 
though this nucleus is naturally stable (negative Q-value) against only light clusters with 
masses A<i < 12, the calculated cluster decay half lives for A<i > 12 are also large enough 
(Tim > 10 80 s) to term this nucleus as a stable nucleus against all cluster decays. This kind 
of stability could apparently be due to stable deformed shell closures at Z=N=38, predicted 
earlier in many other calculations [2,3]. Alternatively, if these nuclei are prepared in heavy-ion 
collisions, then, depending on the excitation energy of the compound nucleus formed, both 
fission (also called, fusion-fission) and cluster decay are the viable processes. The present day 
experiments are directed at these studies (see e.g. [4-7] and earlier references there in). 

It now seems accepted that compound systems with A < 42 are characterized by nuclear 
orbiting phenomena (the deep inelastic process), although a considerable amount of yield 
due to fusion-fission could not be ruled here too [8,9]. On the other hand, the systems with 
A ~ 47 — 60 are strongly the cases of fusion-fission process (fully energy-damped fragments) 
since for all the cases studied so-far the observed yields are independent of nuclei in the 



entrance channel and no strong peaking of yields is observed near the target and projectile 
masses [4-7]. In all these cases, asymmetric mass splitting is favoured and hence lie far 
below the Businaro-Gallone transition point [10] (the fissility parameter x(= Z 2 /50A) is less 
than the xbg = 0.396 for t = and this value decreases as £- value increases). For nuclei 
with A ~ 80, the situation is not so clear. Only three experiments are made [11,12] that 
form the compound systems 3gSr, to-ZV and f^Kr and one calculation is available for 80 Zr 
[13]. Notice that fissility x =0.370 and 0.312 (< xbg), respectively, for 78 Sr and S3 Kr but 
x =0.40(> xbg) for m Zr. Also, at least 78 Sr and 80 Zr are superdeformed in their respective 
ground-states. The interesting result is: whereas the measured mass spectra of 83 Kr [11] 
and 80 Zr [12] are clearly asymmetric and symmetric, as expected, respectively, but that of 
78 Sr [12] is more asymmetric than symmetric. The last experiment on 28 Si + 50 Cr forming 
the compound system 78 Sr* is made only at one energy (Ei a b = \50MeV) and at one angle 
{^lab = 30°). In view of this result we have chosen 76 SV nucleus with x(=0.38) as well 
as deformation (3 larger than that for 78 Sr (estimated = 0.41 for 78 Sr). Such a study 
has not yet been taken up either experimentally or theoretically. Our calculations in the 
following show a clear asymmetric mass distribution, as expected from the point of view of 
Businaro-Gallone transition [10]. This result suggests that a properly angle intergrated mass 
distribution for 78 Sr (not yet measured) should also be clearly asymmetric since x for 78 Sr 
is smaller than that for 76 Sr. 

Theoretically, both the quantum mechanical fragmentation process [13-16] as well as the 
statistical models [8,11,17-19] have been used to explain the measured mass distributions in 
these reactions. The statistical or compound nucleus model calculations, assuming fusion- 
fission process, are made for three possible cases of (i) two spheres separated by a fixed 
distance d=2fm [11], (ii) saddle point shapes, called transition state model [8], and (iii) 
scission point shapes with d taken as a variable [17-19]. In this last case, the Hauser-Feshbach 
formalism [20] is extended to include the decay fragments heavier than a-particle. Perhaps, 
the preformation factor for different fragments should also be added to this extended Hauser- 
Feshback method (EHFM). In a statistical model all open decay channels are taken to be 



equally populated. This is true as long as one is talking of 7-decay and the light particle 
evaporation of n, p and a-particle. However, once the heavy fragments are also included, 
a-particle is known to compete with some heavy fragments (the exotic cluster decays) and 
hence, a preformation probability factor between a- and the other heavy-fragments should 
come in. 

In this paper, we use the quantum mechanical fragmentation theory, the QMFT [13-16], 
and the cluster decay model [1,21,22] based on QMFT. According to the QMFT, the binary 
fragmentation is a collective mass transfer process where both the light and heavy fragments 
(including the light particles) are produced with different quantum mechanical probabilities. 
Applications of both the fragmenation and cluster decay processes are made in only a few 
cases for the light systems [13-16] and infact it was on the basis of this theory that fusion- 
fission was first proposed by Gupta and collaborators in 1984 as the possible explanation of 
the observed data on fragmentation of light systems. 

In the QMFT, a dynamical collective coordinate of mass (and charge) asymmetry r/ = 
(and r/z = z\+zl ) * s introduced whose limiting values are and 1 i.e. < r/ < 1 [14,23-27]. 
Since the potentials V(R,r]) and V(R,r]z), calculated within the Strutinsky renormalization 
procedure (V = Vldm + 5U) by using the appropriate liquid drop model (for Vldm) and the 
asymmetric two-centre shell model (for 5U), are nearly independent of the relative separation 
coordinate R, R can be taken as a time-independent parameter and hence solve the stationary 
(instead of time-dependent) Schrodinger equation in rj: 



The R- value for light nuclei is fixed at the touching configuration of two nuclei [13-16]: 
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In this approximation, the fragmentation potential 



V( V ) = -B 1 {A 1 ,Z 1 )-B 2 {A 2 ,Z 2 



) + E c + V P , 
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where Bi(Ai,Zi) are the experimental binding energies [28], E c = ZiZ2e 2 /R and Vp is the 
additional attraction due to nuclear proximity potential, given by the well known pocket 
formula of Blocki et. al [29]. The charges Zj in (3) are fixed by minimizing V(rfz), defined 
by (3) without Vp, in r/z-coordinate. The rotational energy due to angular momentum (Vi) 
is not added here since its contribution to the structure of yields is shown to be small for 
lighter systems [15]. Of course, V\ should be added for a comparison of the relative yields. 
It may be mentioned here that in a more realistic calculation, the two-centre nuclear shape 
should be used, instead of eq. (2). One can then trace the actual nuclear shapes involved. 
For a more quantitative comparison, perhaps R = R\ + R2 + d(< 2fm) would be a better 
choice because then one is closer to the saddle shape. 

The numerical solution of (1), on proper scaling, gives the fractional mass distribution yields 
for each fragment as 

Y{A{) =| Mv(Ai)) | 2 y/B^j, (4) 

(i= 1 or 2). For the mass parameters B m we use the classical hydrodynamical model of Kroger 
and Scheid [30]. For two touching spheres, this model gives a simple analytical expression, 

B m = ±2S ( !Sfi±« - 1) (6o) 
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with 

"=^( 2 -f-S» <»> 

v c = 7rR 2 c R t , R c = 0.4i? 2 (5c) 

and vt = V1 + V2, the total conserved volume. Also, R2 « R\ and R c {^ 0) is the radius of a 
cylinder of length R t , whose existence allows a homogenous, radial flow of mass between the 
two fragments. Here m is the nucleon mass. 

The nuclear temperature effects in (4) are also included through a Boltzmann-like function 

00 

I 1>R P= E I # I 2 exp(-E^/e) (6) 
with 6, the nuclear temperature in MeV, related to the excitation energy as 

E* = -A6 2 - 6. (7) 
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Furthermore, in some of the calculations here, temperature effects are taken to act also on 
the shell effects as follows [31] 

V = V LDM + SUexp(-e 2 /2.25). (8) 

Similarly, mass parameters should also vary with temperature but no useable prescription 
is available todate. A constant average mass is taken to mean a complete washing of shell 
effects in it. 

For cluster-decay calculations, we use the preformed cluster model (PCM) of Malik and Gupta 
[21,22]. This model, based on the QMFT, also uses the decoupled approximation to R- and 

77-motions and define the decay half-life Ti or the decay constant A as 

2 

A = ^ = Po-P- (9) 
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Here, Po is the cluster preformation probability at a fixed R, given by the solution of the 
stationary Schrodinger equation (1). At R = R\ + R2, Pq = Y(A2), given by eqs. (4) and 
(6). For ground-state to ground-state decay, v = 0. In the following, we choose R = R\ + R2 
(instead of R = Ro, the compound nucleus radius, where V(Ro) = Q-value) since this 
assimilates the effects of both deformations of the two fragments and neck formation between 
them [32]. 

P is the tunnelling probability, which can be obtained by solving the corresponding stationary 
Schro- dinger equation in R. Instead, Malik and Gupta calculated it as the WKB penetrability 
which for the tunneling path shown in figure 1 of [22] is given by 

P = PiPb (10a) 

with 



P i = exp(-l[* {2 f i[V(R)-V(R i )]} 1 / 2 dR (106) 
P b = exp{-- / {2fi[V(R) - Q]} 1/2 dR. (10c) 



This means that tunneling begins at R = R\ + R2 and terminates at R = R\, with V(Rb) = Q- 
The de-excitation probability between Pi and is taken to be unity here. Both (10b) and 



(10c) are solved analytically [21,22]. Apparently, we are considering here the decay from the 
ground-state of the parent nucleus to the ground-states of the decay products. On the other 
hand, if the compound system is excited or the system ends in the excited state of one or 
both the decay products, the Q-value has to be adjusted accordingly (as discussed in the 
following) . 

In eq. (9), v is the assault frequency, given simply as 

■>= ^ = mi^iR», (id 

where Ei = (A\/A)Q is the kinetic energy, taken as the Q-value shared between two frag- 
ments, and n = m( Al ^ 2 ) is the reduced mass. 

Figure 1 shows the fragmentation potential for 76 Sr, plotted as a function of the light fragment 
mass A2. Calculations are made in steps of one nucleon transfer. We notice in Fig. 1 that 
potential energy minima lie only at N=Z, A=4n nuclei, showing the strong shell effects of 
a- nuclei. It is important to realize here that this fragmentation potential is independent of 
the nuclei in the entrance channel and the light (^2) and heavy (Ai = A — A2) fragments 
occur in coincidence. 

The calculated fractional yields are shown in Fig. 2 for the fission of Sr from ground-state 
[y = 0) and at an arbitrary ~ 32 MeV of excitation energy (9 = 2 MeV). We notice that 
in the ground-state, almost all the yield is taken away by the a-particle alone, but as the 
compound system is heated up, other fragments also show up. The interesting results are 
(i) the yields are largest for a-nuclei, and (ii) the mass distribution is strongly asymmetric. 
The fact that both the above mentioned results also hold good in the very small yield zone 
of symmetric fragmentation, we have shown in Fig. 3, the renormalized fractional yields, 
calculated for the fragment masses 15 < Ai < 61 only. Apparently, the a-nuclei are favoured 
and the overall mass distribution is asymmetric. Also, the role of temperature in increasing 
the yields, more so for the symmetric and nearly symmetric fragments, is also shown in Fig. 
3. 

We have also analyzed the role of temperature on shell effects SU. For this purpose, we have 



redone our calculations by using the theoretical binding energies [33] where Vldm and 5U 
contributions are tabulated separately. The calculated fragmentation potentials at different 
temperatures and the resulting yields are shown, respectively, in Figs. 4 and 5. We notice 
that temperature effects are large but the general character of the mass distribution remains 
unaffected. Our earlier calculations [15,16] show that the contribution of shell effects in mass 
parameters B m are also of similar orders and act in the same way i.e. without disturbing the 
general character of the mass distribution. Following [16], the small entrance channel effects 
in the observed mass distributions could be assimilated by the empirically fitted B m (r]). 

Our calculated decay half-lives for the decay of Sr from its ground-state to the ground- 
states of all a-nuclei clusters with positive Q- value are given in Table 1. Of course, the decay 
could also occur into the excited states of the daughter and/or cluster nuclei. This would 
mean decreasing Q-value to Q — Ei, where, say, E; L is the energy of first excited state of 
daughter nucleus. This type of decay is studied in the following paragraph (Table 2). We 
notice in Table 1 that the penetrabilities P are very small. This is particularly so for 12 C 
decay since its Q-value is relatively very small. Apparently, all the predicted decay half-lives 
are very large (Ti/ 2 > 10 80 s) and 76 Sr nucleus can be said to be stable against all possible 
cluster decays. Another interesting point to note is that the daughter nuclei for 20 Ne- and 
32 S-decays are the doubly magic N=Z nuclei. Hence, in view of the so-far observed cluster 
decays [34] , these two decays should be the most probable decays. However, the preformation 
factors Pq for both 20 Ne and 32 S clusters are also very small, though the decrease of Po with 
the increase of cluster size is very much in agreement with the situation for observed decays 
(see e.g. Fig. 13 in the review [34]). 

The role of excitation energy E* of the compound system is presented in Table 2. Once again 
the decays could end into the ground-states or excited states Ei of the fragments. For decay 
into ground-state(s) of the fragment(s), the effective Q-value (Q e ff) will become Q + E* but 
if it goes into an excited state E\ of one fragment, Q e ff = Q + E* — E±. Both the cases 
are studied in Table 2 for some arbitrary values of E* and the first excited states of the 
daughter fragments. Interesting enough, the decay constant A increases (or T X j 2 decreases) 



considerably. Such calculations for light nuclei are made for the first time and have become 
relevant because it is now possible to study experimentally [35] the above mentioned fine 
structure effects of decay into the excited states of fission fragments. Once the data becomes 
available, it will also help deciding between the fission and cluster decays of these light nuclear 
systems. 

Summarizing, we have presented our calculations for fission and cluster decays of 76 Sr. Both 
the cases of 76 Sr in the ground-state and produced as an excited compound system in heavy- 
ion reactions are studied. Here fission is treated as a collective mass transfer process and 
cluster decay studies are based on a model allowing preformation of clusters. Calculations 
show that, in a clear asymmetric mass distribution, w Sr nucleus allows preferential a-nuclei 
transfer resonances as well as decays. The only experimental study available on this system 
is la and 2a transfer products from 36 Ar on 40 Ca at an energy near the Coulomb barrier 
(V c = 53.6 MeV) [36]. Our calculations suggest that for obtaining the complete fission 
products, one has perhaps to go to at least double the Coulomb barrier energies. Also, 
following Sobotka et. al [11], use of inverse kinematics (projectile heavier than the target) 
may be an additional help. This provides a large center-of-mass velocity which facilitates 
the verification of full momentum transfer and easy identification of the fragment's atomic 
number at higher incident energies. Also, the high-enegy solution at forward angles should 
enhance the observation of compound-nucleus decay and virtually eliminate any possible 
deep-inelastic contribution. 
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Figure Captions: 

Fig. 1. The mass fragmentation potential for 76 SV, calculated by using the experi- 
mental binding energies. Only the light fragments A2 are shown along the x-axis. 
The other fragment A\ = A — A2 ■ 

Fig. 2. The calculated mass distribution yields for the fission of 76 Sr from the ground- 
state and at a temperature 6 = 2 MeV (E* « 32 MeV). The range of fragment 
masses is 1 < Ai < 75. 

Fig. 3. The same as for Fig. 2 but at 9 = 2 and 3 MeV and for the range 15 < Ai < 61. 
The fragmentation potential used is the same as in Fig. 1 but for the range 
15 < Ai < 61 only. 

Fig. 4. The same as for Fig. 1 but by using the theoretical binding energies and for the 
mass range 15 < Ai < 61 only. Also, the temperature effects on shell corrections 
5U are shown for 9 = 2 and 3 MeV. 

Fig. 5. The same as for Fig. 3 but by using the potentials of fig. 4. 



